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Abstract. The large elliptic flow observed in Au-Au collisions at RHIC is often 
put forward as a compelling evidence for the formation of a strongly-interacting 
, quark-gluon plasma. The main argument is that the measured elliptic flow 

I> ' is as large as the value given by fluid-dynamics models that assume complete 

04 , thermalization. It is argued that this claim may not be justified, since a detailed 

examination of experimental data rather suggests that the system created is 
. not fully equilibrated at the time when anisotropic flow develops. 
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q ' 1. Introduction 

£h ■ One of the salient results of the heavy-ion programme at RHIC is the measurement 

with unparalleled detail and accuracy of the anisotropy in the transverse-momentum 
distributions of particles emitted in the collisions, the so-called "anisotropic flow" . 
r> | The four Collaborations have provided plenty of data on the first harmonics v\ ( "di- 

rected flow"), V2 ("elliptic flow") and V4 in the Fourier expansion of the azimuthal 
distribution of particles, as a function of the particle transverse momentum and 
rapidity, for various particle species [|TJ |21 El El These data triggered immense 
interest, as it was claimed that for the first time they could be reproduced, together 
with the transverse-momentum spectra of identified particles, by hydrodynamical 
models Such an agreement supposedly necessitates that the matter created in 
the collisions be thermalized after about 0.6 fm/c, and that its equation of state be 
soft. 

The claim led to a huge amount of theoretical studies which investigate whether 
such a short thermalization time can be accounted for in microscopic models of the 
collision Meanwhile, more phenomenological studies are still needed, to question 
the uniqueness of the interpretation of the experimental findings. It actually turns 
out that there is ample room for an alternative reading of the anisotropic-flow data, 
provided one assumes that equilibration in the collisions is incomplete [|7]. In that 
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view, I shall first recall in Sec. [21 various definite predictions of ideal fluid dynamics 
regarding anisotropic flow. The contrasting predictions of an out-of-equilibrium 
scenario will then be presented in Sec. |3 in particular, it will be shown that the 
latter assumption elucidates in a natural way several features of the data that were 
left aside by the ideal-fluid explanation. 

2. Anisotropic flow in ideal fluid dynamics 

In this Section, I shall list different firm predictions for anisotropic flow derived 
within ideal fluid dynamics. Before that, let me first briefly recall the physical 
prerequisites for using a hydrodynamical description of the evolving matter in heavy- 
ion collisions. 

2. 1 . Ideal fluid dynamics: the basic physics ingredients 

In fluid dynamics, it is customary to sort fluid flows into different categories ac- 
cording to their physical properties, using dimensionless numbers. Thus, viscous 
(resp. inviscid, also referred to as "ideal") flows are characterized by small (resp. 
large) Reynolds numbers Re = eLv^n/rj, where s, Wfl u id and r\ are the energy 
density, velocity and shear viscosity of the fluid, and L some characteristic length 
in the system. Similarly, the Mach number Ma = vft U id/c s , where c s is the speed of 
sound in the fluid, quantifies the difference between incompressible (Ma <C 1) and 
compressible (Ma ~ 1) flows. Finally, the Knudsen number Kn = X/L — where 
A is a mean free path — marks the disparity between systems with low numbers 
of collisions per particle (large Kn), which behave like free-streaming gases, and 
liquid- like systems in which each particle experiences many collisions (Kn < 1). 

The three above-mentioned numbers are actually related to each other: since 
1] — ec s X, one finds at once Ma = Kn x Re. Now, the matter created in heavy-ion 
collisions expands into the vacuum, hence the corresponding flow is compressible 
and Ma is of order unity. This implies that the expansion of the created fluid 
satisfies the relation Re ~ 1/Kn: if one can show that Kn is small, it means that 
the fluid viscosity is small, in accordance with the "ideal liquid" paradigm 

In the following, I shall present results on anisotropic flow in both regimes where 
Kn ~ 1 fSec. and Kn <§; 1 (Sec. 12. 2|) . In the latter case, which will be referred 
to as the "ideal fluid" , the large number of collisions per particle leads to some local 
thermal equilibrium. One distinguishes several types of equilibria: cither "kinetic" , 
i.e., equilibrium with respect to elastic collisions, or "chemical", with respect to 
inelastic interactions. It is important to keep in mind that these two equilibria 
do not necessarily hold simultaneously. They are also probed by different observ- 
ables, namely the relative abundances of particle species for chemical equilibrium, 
while the constrains on momentum distributions imposed by kinetic equilibrium 
are rather investigated with the help of anisotropic flow and HBT correlations 
Unless stated explicitly otherwise, any reference to equilibrium or equilibration in 
the remainder of this paper will actually only concern kinetic equilibrium. 
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Impact-parameter dependence of the elliptic-flow buildup 
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Fig. 1. Time evolution of t^/e for different values of the impact parameter, for a 
gas of massless particles (c s = 1/V3). 

2.2. Predictions for anisotropic Sow 

Following the original predictions in Ref. [ , the use of hydrodynamics provides 
a simple intuitive picture for the physics of anisotropic flow. The initial spatial 
anisotropy in the transverse plane of the overlap zone between two nuclei in a non- 
central collision results in a stronger pressure gradient in the direction of impact 
parameter ( "in-plane" ) than perpendicular to that direction ( "out-of-plane" ) . As a 
consequence, in-plane particles acquire more momentum than out-of-plane particles, 
leading to an anisotropy of the transverse- momentum distributions. 

More quantitative statements can be made, based either on analytical calcula- 
tions or on Monte-Carlo computations. Thus, simulations of the development of the 
average elliptic flow v 2 show the existence of various scalings [|7]- As a first example, 
the time development of v 2 is independent of the centrality of the collision: if one 
scales v 2 by the initial spatial eccentricity e = (y 2 — x 2 )/(y 2 + x 2 ) and studies how 
it evolves with time measured in units of R/c s , where R — (l/(x 2 ) + l/(j/ 2 )) -1 ' 2 
quantifies the size of the overlap region, one finds a universal curve for most values 
of the impact parameter, except for the most peripheral collisions (see Fig.EJ. This 
universal behaviour shows in particular that the typical buildup time for elliptic 
flow is R/c s , i.e., about 2-4 fm/c for Au-Au or Cu-Cu collisions, in agreement with 
the findings of transport-model computations [ 111! I12| . 

While Fig. ^ shows that the final v 2 value depends on the shape of the over- 
lap region, v 2 oc e, the scale invariance of ideal fluid dynamics implies that it is 
independent of the system size, i.e., of R. Note, however, that this system-size 
invariance does not allow a straightforward extrapolation from one system to the 
other (say from Au-Au to Cu-Cu collisions), because the initial conditions do not 
scale accordingly (different nuclei have different density profiles). 

For given system size and shape, the final elliptic-flow value varies with the 
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speed of sound c s . More quantitatively, if one assumes a constant c s throughout 
the system evolution, then the final t>2 increases with c s as soon as c s > 0.1 — it 
even becomes proportional to c s for c s > 0.3 [\7\. 

To close this Section, let me mention a few results that were obtained analyti- 
cally, exploiting the fact that the ideal-fluid assumption is equivalent to considering 
the limit of small freeze-out temperature. In Ref. [ I13| . it was emphasized that 
emitted particles fall in a natural way into two categories, namely "slow" particles, 
defined as those whose velocity equals that of the fluid at some point on the freeze- 
out hypersurface, and "fast" particles, which are faster than the fluid at freeze-out. 
For both categories of particles, definite predictions regarding anisotropic flow can 
be made. Thus, the dependence of elliptic flow on the particle velocity «2 (pt/w) > 
where pt and m are the particle transverse momentum and mass, should be identi- 
cal for all types of slow particles (except pions, whose mass is not much larger than 
the freeze-out temperature). The same property holds for all other flow harmonics 
v n (pr/m). This in particular implies a mass-ordering of ^(pr), the heavier par- 
ticles having smaller flow at a fixed transverse momentum. It turns out that the 
mass-ordering of V2{pt) also holds for fast particles. For the latter, it was shown 
that the various even flow harmonics are related, the most important relation being 
va(pt) = i"2(pt) 2 /2, valid for each particle species in each rapidity window provided 
the transverse- velocity profile at freeze-out is not too different from an ellipse [I13j. 

3. Out-of-equilibrium scenario 

Let me now turn to the case in which the mean number of collisions per particle, 
Kn^ 1 , is insufficient to lead to any (local) equilibrium. 

3.1. Anisotropic How in the out-of-equilibrium regime 

Determining the precise dependence of anisotropic flow on Kn~ x (as well as on the 
collision cross sections) requires a transport model which is beyond the scope of the 
present study, yet general predictions are nonetheless feasible [[?]• 

Thus, it is natural to expect that elliptic flow increase with the number of col- 
lisions, since V2 vanishes in the absence of final-state interaction, whereas in the 
l&rge-Kn , ideal-fluid limit V2 is finite. This growth, however, eventually satu- 
rates for some Kn~ x value, for which the system equilibrates (see Fig. [3 ideal-fluid 
dynamics is expected to yield the maximum V2, as hinted at by transport compu- 
tations, which always give smaller values [ I14j). The value at which V2 saturates 
obviously depends on the system shape, i.e., on the initial eccentricity e, since i>2 oc e 
in hydrodynamics; however, the corresponding value of the Knudsen number should 
be quite independent of e, as individual parton-parton or hadron-hadron interac- 
tions do not know anything about the geometry of the nucleus-nucleus collision. 
As a consequence, the slope in the region where V2(Kn ) grows should be roughly 
proportional to e; conversely, the ratio t>2/e should be (almost) independent of the 
centrality of the collision for a fixed value of Kn~ x . 
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Fig. 2. Sketchy representation of the variation of V2 with the mean number of 
collisions per particle Kn^ 1 . 

The same reasoning applies to U4, which should also increase with the mean 
number of collisions, and then saturate (at a value oc e 2 ). If V4 does not grow much 
faster than u 2 , then the ratio v^jv\ is a decreasing function of Kn^ 1 , which reaches 
a minimum in the ideal- fluid limit. Since the hydrodynamical prediction for the 
ratio is 1/2, then V4 > v%/2 when equilibrium is not reached [\7\. 

In the out-of-equilibrium regime, both vi and W4 increase with the number of 
collisions. Now, Kn = L/X may vary for two reasons, due to changes either in 
the system size L, for which a natural choice is L — R, or in the mean free path 
A. The former possibility strongly contrasts with the scale invariance of anisotropic 
flow within the ideal- fluid description [ Ej . In turn, the mean free path A = 1/na 
depends on both particle density n and cross section cr, and any change in one of 
these will affect the Knudsen number, hence the flow coefficients when the system 
is not equilibrated. 

To summarize this part, in an out-of-equilibrium scenario the anisotropic flow 
coefficients vi , 1*4 depend significantly on the system size R (even if the system shape 
e is fixed), on the particle density and on the interaction cross section, and they are 
related by v^jv\ > 1/2, the ratio increasing as one goes further from equilibrium. 

3.2. Confronting RHIC data and the out-of-equilibrium scenario 

Let me now show that available flow data support the idea that the system created in 
Au-Au collisions at RHIC is not equilibrated, at least at the time when anisotropic 
flow develops. 

The first element that supports the out-of-equilibrium scenario is its ability to 
explain the rapidity dependence of elliptic flow [E3 ^j. Thus, the fact that V2(y) 
follows closely the rapidity distribution dN / dy from midrapidity up to the forward 
regions, where both exhibit the "limiting fragmentation" property across different 
beam energies [El, is naturally explained within a model where v-i varies with the 
number of collisions per particle. Now, the identity (cR/c s ) n(R/c s ) = (1/5) dN/dy 
yields the particle density at the time when anisotropic flow builds up, which gives 
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where 5 measures the transverse area of the collision zone. Equation shows that 
the incomplete-equilibration prediction v-i oc Kn~ l translates into 1*2(2/) °c dN/dy. 
On the other hand, the few 3-dimcnsional hydrodynamical models cannot reproduce 
the data, either overestimating 1*2(2/) a t intermediate rapidities y ~ 3 [ I15j or, when 
those values are well predicted, underestimating the midrapidity elliptic flow [I17|. 

A similar instance of experimental result which finds a natural explanation in 
a non-equilibrium model is the growing discrepancy between V2{pt) and the ideal- 
fluid prediction as transverse momentum increases [I18|. As a matter of fact, with 
increasing px the interaction cross section is expected to decrease, diminishing the 
mean number of collisions per particle Kn" 1 oc cr, which is thus increasingly further 
from the number needed to ensure equilibration. In turn, this leads to an increase 
of the difference between the hydrodynamical t>2 and the out-of-equilibrium value. 

One could claim that the previous two arguments concern only a negligible frac- 
tion of the particles, while the bulk, at midrapidity and moderately low transverse 
momenta, would still be in equilibrium. If this were true, then the elliptic flow 
V2 at midrapidity, averaged over transverse momentum, should be proportional to 
the eccentricity. However, it turns out that the ratio i^/e is not constant across 
ccntralitics in Au-Au collisions the data do not exhibit the scale invariance of 
ideal fluid dynamics. On the contrary, the values of vijt rather seem to be scaling 
linearly with the control parameter (l/S) dN/dy [EH> i- e -i according to Eq. and 
assuming that c s and the cross section remain roughly constant for the various cen- 
tralities, with the number of collisions. This proportionality, without a single hint 
at any saturation, even extends down to the values measured in Pb-Pb collisions 
at the CERN SPS, lending further credence to the out-of-equilibrium scenario. 

Yet another indication that the system created in Au-Au collisions at RHIC 
is not equilibrated is provided by the ratio v^jv\. STAR and PHENIX reported 
values which are significantly larger than the ideal-fluid value of 1/2 [El I2()| . 

As a final evidence that approaches based on equilibrium assumptions may not 
be appropriate at RHIC, let me point out another failure of the ideal- fluid descrip- 
tion, which can be resolved provided one drops some equilibrium constraint. Thus, 
although hydrodynamics describes properly the transverse-momentum dependence 
of v% of identified particles in minimum bias collisions, on the other hand it misses 
the centrality dependence of 1*2 (pt): the ideal-fluid calculation that reproduce the 
minimum-bias average overestimate elliptic flow in peripheral events and underesti- 
mate its value for pions in central collisions [E] . Even though the former discrepancy 
is not really surprising — one does not expect the hydrodynamical description to 
hold for peripheral collisions, where the system size is too small to allow any equi- 
libration — , however the fact that data largely overshoot the so-called ideal-fluid 
limit in central events is a much more serious issue. A plausible explanation could 
be that the analysis of the data for these centralities is not fully reliable, which is 
possible since this is where the anisotropic-flow signal is smallest, hence most diffi- 
cult to extract accurately. Another possibility is that the central data could indeed 
be described by ideal fluid dynamics, albeit with a stiffer equation of state, i.e., a 
larger speed of sound, than what is currently used. This can be done, provided 
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one realizes that the supposedly strong constraint on c s arising from fits to rapidity 
distributions exists if, and only if, one assumes that the system has reached not only 
kinetic equilibrium, but also chemical equilibrium. Once this assumption (which is 
only supported by the success of fits from "thermal" models to particle- abundance 
ratios, as seen also in e + e~ collisions) is dropped, the one-to-one relationship be- 
tween particle number density, which gives the particle distribution, and energy per 
particle, which is responsible for flow, disappears, and the constraint on c s is lifted. 
Allowing now kinetic equilibrium only, one may describe central collisions with a 
hydrodynamical model — however, even if it proved valid for central events, the 
ideal-fluid description would not hold in more peripheral collisions. 

It could be hoped that a discriminating test between the early-thermalization 
and out-of-equilibrium scenarios would be provided by measurements of anisotropic 
flow in Cu-Cu collisions [u\: the change in system size would probe the scale in- 
variance of ideal fluid dynamics, which is broken in the absence of equilibrium. 
Unfortunately, the preliminary 1)2 results presented at Quark Matter '05 were quite 
inconclusive, as the values of the three different experimental Collaborations were 
not compatible with each other. These discrepancies call for further investigations 
of the possible sources of systematic error on the measurements, in particular fluctu- 
ations of the signal (whether they arise from truly physical effects, or from binning 
issues) and non-flow effects. The latter were shown to affect the flow analysis in 
Au-Au collisions, and should be even more important in the smaller Cu-Cu system; 
several new methods of measurement that are free from their bias were specifically 
introduced [ I21II22| and could be used. Unfortunately, disentangling non-flow effects 
from fluctuations of the flow itself is not a trivial task [I23|. 

Measurements of anisotropic flow with great accuracy will also be needed to 
investigate another expected behaviour, namely that the ratio v^jv\ increases at 
high px and/or when going away from rapidity, since in both these regimes the 
measured trends of V2(pr) and V2(y) suggest that Kn~ x is decreasing. 

Eventually, one can anticipate that anisotropic flow measurements in Pb-Pb 
collisions at ^[sWn = 5.5 TeV at the CERN LHC will help confirm the out- 
of-equilibrium interpretation of RHIC data. Thus, unless the ideal-fluid limit is 
marginally reached in the most energetic central Au-Au events at RHIC, the aver- 
age V2 at midrapidity, scaled by the spatial eccentricity, should be larger at LHC 
than it is at RHIC. Conversely, the ratio V4(pt)/v2(pt) 2 should be smaller, ap- 
proaching the hydrodynamical value of 1/2. 

4. Conclusion 

In summary, I have shown that dropping the assumption of kinetic equilibrium at 
the time when anisotropic flow develops allows me to describe in a satisfactory, 
consistent manner several features in the flow data that cannot be accommodated 
in ideal-fluid models. The only instance where a hydrodynamical approach may 
remain appropriate is in central events; however, even though kinetic equilibrium 
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could be reached, the constraint of chemical equilibrium has to be abandoned if one 
is to reproduce the i>2 data. 

Although they permit a better description of the data, the predictions of the 
out-of-equilibrium scenario presented here are admittedly quite crude, and deserve 
further dedicated studies. A transport model would allow a more quantitative inves- 
tigation of how the system evolves from an non-equilibrated state to a thermalized 
one, answering various questions as when (for which value of the Knudsen number) 
does anisotropic flow reach the hydrodynamical limit? How far are RHIC Au-Au 
data from this limit? In that respect, the ratio v^jv\ might be a more sensitive in- 
dicator than V2 itself, as it seems to be further away from the ideal-fluid prediction. 
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